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SUMMARY 


Mat hen at icai  relationship  for  the  dielectric  and  permeability  constant  of 
single  particle  and  composite  materials  have  been  formulated  and  tested  in 
experimental  conditions.  For  dilute  concentrations  of  particles  imbedded  in 
cinder  materials  reasonable  agreement  was  found  between  measured  and 
calculated  values  of  this  constants.  The  theoretical  formulation  is  based 
upon  a  two  parameter  theory  (magnetization  and  magnetic  anisotropy  field). 


EFFECTS  OF  DEPOLARIZATION  ON  THE  MICROWAVE  PROPERTIES  OF 

COMPOSITES. 


C.  Vittoria 

Northeastern  University 
Boston,  MA  02115 


1  BACKGROUND 

The  microwave  properties  of  composite  materials  are  of  great 
interest  to  many,  since  it  is  possible  to  control  or  design  the 
constituents  of  a  composite  in  order  to  obtain  a  desireable 
property  of  the  composite.  There  are  many  physical  parameters 
that  can  effect  the  microwave  properties  of  a  composite  and  they 
are:  binder  material,  particle  size,  shape,  electrical  and 
magnetic  properties  of  the  particle  and  binder,  density  of 
particles,  randomness  of  particle  distribution,  etc...  In  this 
report  we  will  concentrate  our  attention  on  the  effect  of 
particle  shape  on  the  microwave  properties  of  a  composite. 

In  particular  we  will  consider  metallic  as  well  as  insulating 
particles  distributed  randomly  in  a  known  binder  material.  The 
effect  of  the  particle  shape  in  composite  is  to  induce  electrical 
and  magnetic  depolarization  fields  within  the  particle  itself  as 
well  as  beyond  the  particle.  Depolarization  fields  are  induced 
within  magnetic  particles,  if  the  magnetization  is  uniform 
throughout  the  particle  -  single  domain.  In  multi-domain 
configurations  the  depolarizing  field  within  the  magnetic 
particle  is  reduced  and  may  be  assumed  to  be  zero  for  practical 
purposes.  The  density  of  magnetic  domains  nucleated  in  a 
particle  depends  on  the  size  and  shape  of  the  particle.  Similar 
observations  can  be  made  for  the  electric  depolarization  field, 
if  the  particle  is  ferro-electric,  for  example.  However,  for 
most  applications  scalar  dielectric  constant  materials  are  used 
in  composites.  Hence,  we  will  assume  scalar  dielectrics  in  this 
report.  In  such  cases  the  shape  of  the  particle  plays  a  major 
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role  in  the  magnitude  of  the  depolarizing  electric  field  within 
the  particle.  For  dilute  concentrations  of  particles  in  a  binder 
it  is  meaningful  to  confine  depolarizing  fields  within  the 
particle.  This  assumption  breaks  down,  if  the  density  of 
particles  increases  beyond  25%  (by  volume)  of  the  composite 
material.  We  will  now  consider  firstly  metallic  iron  alloy 
particles  with  dimensions  of  the  order  of  20^®  imbedded  in 
castolite.  For  the  insulating  particles  we  will  consider  ferrite 
particles  imbedded  in  a  binder  material  characterized  by  high 
dielectric  constant. 


2  METALLIC  PARTICLES 


Composite  materials  were  fabricated  in  which  iron  alloy  particles 
were  imbedded  into  castolite  binder.  Castolite  is  a  dielectric 
material  with  a  corresponding  value  of  dielectric  constant  of 
2. 59- jO.  05  over  the  frequency  range  of  1  to  20  Ghz.  The  iron 
alloy  particles  were  sifted  through  a  20Mm  grid  mesh.  Thus,  some 
particles  were  greater  and  others  less  than  20»m  in  size.  The 
shape  of  the  particles  was  unknown.  The  samples  were  designated 
as  one  set  being  greater  and  the  other  less  than  20^m.  Toroids 
of  90%  castolite  and  10%  particles  (by  volume)  were  fabricated  in 
the  microwave  laboratory  for  the  purpose  of  measuring  t 
(dielectric  constant)  and  n  (permeability  constant)  of  the 
composite.  A  network  analyzer  was  used  to  measure  e  and  ^  from 
0.5  to  18.6  GHz. 

For  the  composite  (characterized  by  particles  greater  than  20^m 
in  size)  the  average  dielectric  constant  was  measured  to  be 

€  av  =  5.40  -  jO. 080  ( >2 0^m)  , 


whereas  for  the  other  case 

=  4.25  —  jO. 075  (<20^m) 
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The  value  of  «  remained  fairly  constant  from  0.5  to  18.6  GHz. 

In  analyzing  the  data  we  made  use  of  the  well  known  formula  of 
Maxweli-Garnett ( 1  ) 


,  F  (e  -  Q _ 

fav  eb  (1  -  F) ( 1  +  L( t  -  «b))  +  F 


(1) 


where  t  =  dielectric  constant  of  particle 

=  dielectric  constant  of  binder  material 
F  =  fractional  volume  of  particle  loading 
L  =  depolarizing  factor 


For  a  conductive  medium  «  may  be  related  to  the  conductivity  by 
the  relation 

«  =  a/ jw 

a 

o 

where  o  —  - 

1  -  j  ur 


r  -  relaxation  time  between  collisions. 
a 0  *  DC  conductivity. 


If  the  particle  is  conductive,  eg.  (1)  reduces  to  approximately 


e 

av 


“aTL 


+  f  b 


(2) 
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Re-write  above  equation  as  follows,  since  F  =  0.10  for  our 
composites . 


fb 


-  1 


(3) 


Figure  (1)  shows  a  plot  of 
our  case  that  2  =  € ^ 


fav/£b  versus  ‘  We  see  that  for 
for  the  measured  value  of  ‘avAir,  a  2 


corresponding  the  composite  containing  particles  greater  than  2^m 
and  2  =  for  the  other  composite  (<  20nm) . 

The  fact  that  the  depolarizing  factor  depends  on  the  size  of  the 
particles  implies  that  the  shape  of  these  particles  are  dependent 
on  the  size  of  the  particles.  The  particles  tend  to  be  more  like 
disc  types  for  particles'  diameter  bigger  than  20/iin. 

The  same  depolarizing  factor  must  also  be  effective  in  describing 
the  magnetic  properties  of  the  particle.  Hence,  if  we  were  to 
apply  equation  (1)  to  permeable  materials,  it  would  be  written  as 
follows 


av 


1  + 


(1 


f(m  -  l) _ 

F)  ( l  +  Hu  -  D)  +  F 


(4) 


where  >ib  ■  permeability  constant  of  binder  =  1 
jj.  »  permeability  constant  of  particle. 


The  permeability  of  the  particle  may  be  approximated  as  follows 
for  relative  large  size  samples 


M  ,  i  +  x*Jd  +  JX/2Q) 

1  -  x2  +  jx/Q 


(5) 
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where 


Xo  =  4^M/Ha 

Q  =  f0/Af 
X  =  f/fG 

4*M  =  Saturation  Magnetization. 

Ha  =  Uniaxial  magnetic  anisotropy  field. 
f0  =  ferromagnetic  resonant  frequency  =  7  HA . 

f  =  Operating  frequency. 

Af  =  ferromagnetic  resonance  linewidth  7  2*M 

7  «■  2 . 8MHz/0e  . 

The  factor  Q  is  a  measure  of  the  material  losses  at  microwave 
frequencies.  It  is  important  to  note  that  the  above  equation  is 
only  appropriate  to  particles  in  which  multi-domains  are 
nucleated  in  a  fairly  high  density.  Particles  in  the  order  of 

20pm  contain  a  very  high  density  of  multi-domains.  It  is 
estimated  that  one  needs  to  fabricate  particles  in  the  order  of 
1/4  pm  in  order  to  energetically  favor  single  domain  excitations 
instead  of  multi-domains.  For  single  domain  excitations  equation 
(5)  is  inappropriate,  because  eq.  5  does  not  take  into  account 
of  demagnetizing  fields  or  depolarizing  fields.  There  are 

other  complications  which  the  reader  should  be  aware  of  in  the 
use  of  equation  (5)  .  In  a  conductive  medium  eddy  currents 
effects  tend  to  mask  or  screen  the  external  electromagnetic 
fields.  As  such  there  are  diamagnetic ( 3 ^  contributions  to  p .  In 
addition  the  electromagnetic  wave  is  attenuated  in  the  metallic 
particle.  The  attenuation  is  further  enhanced  by  the  exchange 
conductivity  effects^4).  The  result  of  this  is  to  slightly  shift 
f0  from  the  simple  relation  7HA  and  increase  Af.  Estimates  of 
these  shifts  can  be  calculated ^ 4 )  and  incorporated  in  the 
expression  for  p  in  eq.  (5)  by  slightly  modifying  fQ  and  Af. 

Let's  calculate  p  for  a  single  particle  explicitly  assuming  the 
following 


4 wM  =  21,500  G 
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Ha  =  500  Oe 
if  =  30  GHz 
f0  =  1.5  GHz 

The  values  for  these  parameters  are  assumed  on  the  basis  that  we 
are  dealing  with  iron  rich  alloy  particles  and  the  cubic  magnetic 
anisotropy  field  is  in  the  order  of  500  Oe  and  4wM  is  close 

to  that  of  pure  iron.  The  calculat ional  procedure  is  as  follows: 
The  above  parameters  together  with  an  assumed  frequency  value  are 
substituted  into  eq.  (5)  and  ^  is  calculated  for  the  single 
particle.  In  order  to  calculate  ^av  of  the  composite  the  value 
of  n  of  the  particle  is  substituted  into  eq.  4  assuming  F  =  0.10. 
In  figs.  (2)  and  (3)  ^av  of  the  composite  are  calculated  for 

<  =  |  and  "  ,  respectively. 

,  ,  / 

The  most  prominent  feature  of  the  calculations  is  that  m  av  is 
constant,  where  ^'av  is  the  real  part  of  Mavr  which  is  in 
contrast  to  measured  values  of  M'av-  In  fig.  (4)  we  have  plotted 
Pav  calculated  from 


'“av  "  1  +  F  (m  -  1)  (6) 

In  equation  (6)  we  make  the  assumption  that  there  are  no 

depolarizing  factors  in  multi-domain  excitations  in  the  magnetic 

particles.  In  fig.  (4)  both  calculated  and  measured  M/av  scale 

inversely  with  frequency  with  approximately  the  same  slope. 

Qquantitatively ,  it  appears  that  inserting  depolarizing  factors 

/ 

in  the  expression  for  n  av  improves  the  fit  to  the  experimental 

observation.  However,  the  agreement  is  only  good  at  low 

/ 

frequencies  and,  furthermore,  calcualted  n  av  is  constant  with 
frequency.  Omitting  depolarizing  factors  eq.  (6)  predicts  a 
scaling  of  ^  av  with  frequency,  but  quantitatively  a  factor  of 
1.7  smaller.  We  believe  that  it  is  important  to  predict  the 

f  , 

frequency  dependence  of  n  av,  since  there  are  other  contributions 
which  can  account  for  the  quantitative  difference.  We  will 
discuss  these  contributions  at  the  end  of  this  section. 
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Also  calculated  M"av  is  in  reasonable  agreement  with  measured 
M"av  considering  the  measured  values  of  M"av  are  relatively  small 
compared  to  one.  Network  analyzer  techniques  are  not  amenable  to 
measuring  small  values  of  ^"av. 

The  difference  between  the  measured  and  calculated  values  of  M"av 
may  be  accounted  by  three  equally  important  factors:  (1)  Since 
the  particles  are  conductive,  only  material  within  the  skin  depth 
distance  from  the  particle  surface  is  magnetically  active. 
Effectively,  this  means  that  the  actual  filling  factor  F  is 
smaller  than  0.1  we  have  assumed  above  in  equation  (6)  toroids. 
(2)  Thus,  if  we  were  to  choose  for  example  F  =  0.04,  better 
agreement  can  be  obtained  between  measured  and  calculated  values 
of  However,  even  assuming  F  =  0.04  there  is  still 
discrepancy  in  M"av-  The  measured  M"av  is  a  factor  of  two  bigger 
than  the  calculated  value.  This  difference  may  be  accounted  for 
by  additional  eddy  current  losses  found  in  conductive  materials. 
This  extra  loss  manifests  itself  in  the  measurement  of  M"av-  (3) 
Finally,  equation  (5)  is  indeed  an  approximate  equation. 
Accurate  expression  for  m  for  specific  multi-domain  conf iguations 
may  be  obtained  in  reference  (6).  The  parameter  values  for  Af 
and  fG  were  chosen  in  a  manner  to  make  the  resonance  curve  highly 
asymmetrical.  Equation  (5)  applies  only  to  resonance  curves 
which  are  symmetrical  with  respect  to  frequencies  below  above  f0 . 
Nevertheless,  we  feel  that  reasonable  agreement  between  measured 
and  calculated  values  of  was  obtained. 

3  INSULATING  PARTICLES 

For  insulating  particles  one  may  ignore  eddy  current ( 3 )  and 
exchange  conductivity^  effects.  With  these  omissions  eq.  (5) 
should  be  more  relevant  to  insulating  particles.  However,  the 
expression  for  e  needs  to  be  modified  to  the  following 
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€  -  £c  1  +  ‘ij2p/((w20  “  w2]  +  3  “  j 


c0  =  dielectric  constant  of  free  space 

cjp  =  plasma  frequency 

ai0  =  electronic  resonance 

Data  was  provided  to  us  on  a  composite  material  which  consisted 
of  the  following  constituents:  Ni-zn  ferrite  particles^7)  which 
exhibited  saturation  magnetization  of  (4*M)  4200  G,  H  =  3600  Oe, 
F  =  0.62,  e  =  (20.9  -  0.48  f)  -  j(0.7  -  0.005f).  The  dielectric 
constant  of  the  binder,  £b,  was  2.8.  In  fig.  (5)  calculated  and 
measured  values  of  £  av  have  been  plotted  as  a  function  of 
frequency  from  2  to  18  GHz.  Calculated  values  were  obtained 
using  equation  (1)  and  assuming  the  following  parameter  values: 
£b  =  2.8,  F  =  0.62  and 


We  now  apply  equation  (5)  and  (6)  to  calculate  ^av  of  the 
composite  assuming  above  magnetic  parameters.  In  fig.  (6)  the 
permeability  of  the  single  particle  is  plotted  as  a  function  of 
frequency.  The  comparison  between  calculated  and  measured  n  is 
reasonable  considering  that  only  two  adjustable  parameters  are 
used  in  eq.  (5)  .  Improved  comparison  may  be  obtained,  if  Af  is 
lowered  with  respect  to  the  assumed  value  of  5.88  GHz.  However 
reduction  of  Af  is  not  consistent  with  maintaining  4*M  constant, 
since  Af  =  72*M.  In  fig.  (7)  we  calculated  ^av  using  eq.  (6) 
together  with  eq.  (5) .  Again  only  two  adjustable  parameters  are 
assumed,  HA  and  4wM,  which  also  give  f0  =  7HA  =  10.08  GHz  and  Af 
-  2jtM  =  5.88  GHz.  In  fig.  7  the  agreement  between  calculated  and 
measured  Mav  is  reasonable  if  we  were  to  lower  fQ  by  as  much  as 
2.2  GHz  in  eq.  5.  There  is  no  physical  reason  to  do  so.  Further 
examination  of  the  material  preparation  of  the  ferrites  it  has 
been  revealed  that  the  ferrites  were  sintered  at  high  temperature 
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prior  to  the  fabrication  of  the  composite.  It  is  well  known  that 
sintering ( 2 ' 5 )  processes  in  ferrites  can  significantly  alter  the 
values  of  HA,  for  example.  Hence,  we  conclude  that  the  value  of 
HA  was  lowered  by  roughly  750  Oe  by  the  process  of  sintering. 

4  CONCLUSIONS 

1.  Depolarizing  factors  are  important  in  the  calculation  of  the 
average  value  of  the  dielectric  constant  in  composite 
materials  consisting  of  dielectric  particles.  For  metallic 
particles  propagation  effects  must  be  considered  especially 
for  composites  made  up  of  relative  large  size  particles  in 
comparison  to  the  skin  depth. 

2.  For  large  size  magnetic  particles  imbedded  in  composites  the 
depolarizing  factor  may  be  omitted  in  the  calculation  of 
the  average  value  of  the  permeability  of  composite 
materials. The  depolarizing  factor  may  be  omitted,  since 
multi-domains  are  nucleated  in  large  size  magnetic 
particles.  The  depolarizing  field  averages  out  to  zero  for 
multi-domain  configuration:  Our  equation  (5)  is  a 
reasonable  estimate  of  n  in  the  regime  of  multi-domain 
excitations.  However,  caution  must  be  exercised  in  the 
use  of  equation  (5),  if  the  domain  density  is  small  and 
close  to  unity.  Exact  calculations  must  be  invoked,  if 
improved  estimates  of  n  are  desireable  for  low  domain 
densities.  For  single  domain  excitations  as  in  very  small 
magnetic  particles  the  reader  is  referred  to  reference  (2) 
for  the  calcualtion  of  /j.  In  this  limit  one  then  must  also 
include  depolarizing  factors  in  the  calculation  of  m  and 
*xav  (eqs.  (5)  and  (4)). 

3.  The  above  two  conclusions  are  applicable  for  dilute 
concentration  of  particles  in  a  composite  material.  For 
large  concentrations  of  particles  in  a  composite,  particle 
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to  particle  interactions  must  be  included.  The  interactions 
may  be  electrostatic  or  magnetostatic  in  nature.  This  is 
still  an  outstanding  theoretical  problem  which  is  of  current 
interest. 

We  have  proposed  a  two  parameter  model  in  the  calculation  of 
both  m  and  Mav,  namely  and  4wM.  f0  and  Af  scale  as  HA 
and  M,  respectively.  However,  microwave  data  on  ^  and  Mav 
seems  to  suggest  at  least  one  more  parameter  is  needed  to 
explain  the  data  completely.  A  simple  way  to  modify  our 
present  analysis  is  to  allow  Af  to  be  a  variable,  when 
substituted  in  eq.  (5).  As  a  future  guide  one  may  assume  a 
maximum  value  of  Af  ^7  2*M  and  a  minimum  value  in  the  order 
of  the  intrinsic  linewidth  of  the  single  particle. 
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Force  Command,  Control,  Communications  and  Intelligence  (C  I)  activities 
for  all  Air  Force  platforms.  It  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineering  support  within 
areas  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other 
ESD  elements  to  perform  effective  acquisition  of  C  I  systems.  In  addition, 
Rome  Laboratory's  technology  supports  other  AFSC  Product  Divisions,  the 
Air  Force  user  community,  and  other  DOD  and  non-DO D  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  re  liability/  maintainability  and  testability. 


